CGG toolkit: Software components for computational genomics.

Public-domain availability for bioinformatics software resources is a key requirement that ensures long-term permanence and methodological reproducibility for research and development across the life sciences. These issues are particularly critical for widely used, efficient, and well-proven methods, especially those developed in research settings that often face funding discontinuities. We re-launch a range of established software components for computational genomics, as legacy version 1.0.1, suitable for sequence matching, masking, searching, clustering and visualization for protein family discovery, annotation and functional characterization on a genome scale. These applications are made available online as open source and include MagicMatch, GeneCAST, support scripts for CoGenT-like sequence collections, GeneRAGE and DifFuse, supported by centrally administered bioinformatics infrastructure funding. The toolkit may also be conceived as a flexible genome comparison software pipeline that supports research in this domain. We illustrate basic use by examples and pictorial representations of the registered tools, which are further described with appropriate documentation files in the corresponding GitHub release.

Targeting oncogenic microRNAs from the miR-371~373 and miR-302/367 clusters in malignant germ cell tumours causes growth inhibition through cell cycle disruption.

BACKGROUND: MiR-371~373 and miR-302/367 cluster over-expression occurs in all malignant germ cell tumours (GCTs), regardless of age (paediatric/adult), site (gonadal/extragonadal), or subtype [seminoma, yolk sac tumour (YST), embryonal carcinoma (EC)]. Six of eight microRNAs from these clusters contain the seed sequence 'AAGUGC', determining mRNA targeting. Here we sought to identify the significance of these observations by targeting these microRNAs functionally. METHODS: We targeted miR-371~373 and/or miR-302/367 clusters in malignant GCT cell lines, using CRISPR-Cas9, gapmer primary miR-302/367 transcript inhibition, and peptide nucleic acid (PNA) or locked nucleic acid (LNA)-DNA inhibition targeting miR-302a-d-3p, and undertook relevant functional assays. RESULTS: MiR-302/367 cluster microRNAs made the largest contribution to AAGUGC seed abundance in malignant GCT cells, regardless of subtype (seminoma/YST/EC). Following the unsuccessful use of CRISPR-Cas9, gapmer, and PNA systems, LNA-DNA-based targeting resulted in growth inhibition in seminoma and YST cells. This was associated with the de-repression of multiple mRNAs targeted by AAGUGC seed-containing microRNAs, with pathway analysis confirming predominant disruption of Rho-GTPase signalling, vesicle organisation/transport, and cell cycle regulation, findings corroborated in clinical samples. Further LNA-DNA inhibitor studies confirmed direct cell cycle effects, with an increase of cells in G0/G1-phase and a decrease in S-phase. CONCLUSION: Targeting of specific miR-371~373 and miR-302/367 microRNAs in malignant GCTs demonstrated their functional significance, with growth inhibition mediated through cell cycle disruption.

RNA sequencing and lipidomics uncovers novel pathomechanisms in recessive X-linked ichthyosis.

Recessive X-linked ichthyosis (RXLI), a genetic disorder caused by deletion or point mutations of the steroid sulfatase (STS) gene, is the second most common form of ichthyosis. It is a disorder of keratinocyte cholesterol sulfate retention and the mechanism of extracutaneous phenotypes such as corneal opacities and attention deficit hyperactivity disorder are poorly understood. To understand the pathomechanisms of RXLI, the transcriptome of differentiated primary keratinocytes with STS knockdown was sequenced. The results were validated in a stable knockdown model of STS, to confirm STS specificity, and in RXLI skin. The results show that there was significantly reduced expression of genes related to epidermal differentiation and lipid metabolism, including ceramide and sphingolipid synthesis. In addition, there was significant downregulation of aldehyde dehydrogenase family members and the oxytocin receptor which have been linked to corneal transparency and behavioural disorders respectively, both of which are extracutaneous phenotypes of RXLI. These data provide a greater understanding of the causative mechanisms of RXLI's cutaneous phenotype, and show that the keratinocyte transcriptome and lipidomics can give novel insights into the phenotype of patients with RXLI.

Small non-coding RNA sequencing reveals global dysregulation of piwi-interacting RNA (piRNA) expression in gonadal malignant germ cell tumours.

BACKGROUND: Analyses of small non-coding RNA (ncRNA) expression in malignant germ cell tumours (GCTs) have focused on microRNAs (miRNAs). As GCTs all arise from primordial germ cells, and piwi-interacting RNAs (piRNAs) have important roles in maintaining germline integrity via transposon silencing, we hypothesised that malignant GCTs are characterised by fundamental piRNA dysregulation. AIMS: We undertook global small ncRNA sequencing in malignant GCTs, in order to describe small ncRNA expression changes for both miRNAs and piRNAs. MATERIALS AND METHODS: We performed small ncRNA next generation sequencing on a representative panel of 47 samples, comprising malignant GCT (n = 31) and control (n = 16) tissues/cell lines. Following quality control and normalisation, filtered count reads were used for differential miRNA and piRNA expression analyses via DESeq2. Predicted mRNA targets for piRNAs were identified and utilised for pathway enrichment analyses. RESULTS: Overall, miRNAs and piRNAs comprised 21.9% and 43.0% of small ncRNA species, respectively. There were 749 differentially expressed miRNAs in malignant GCTs, of which 536 (72%) were over-expressed and 213 (28%) under-expressed. The top-ranking over-expressed miRNAs were exclusively from the miR-371∼373 and miR-302/367 clusters. The most significantly under-expressed miRNAs were miR-100-5p, miR-214-3p, miR-125b-5p and let-7 family members, including miR-202-3p. There were 1,121 differentially expressed piRNAs in malignant GCTs, of which 167 (15%) were over-expressed and 954 (85%) under-expressed. Of note, of the top-20 differentially expressed piRNAs, 16 were over-expressed, of which piR-hsa-2506793 was both top-ranking and most abundant. Mobile element (ME; i.e., transposon)-associated piRNAs comprised 166 (15%) of the 1,121 differentially expressed piRNAs, of which 165 (>99%) were down-regulated. The remaining 955 (85%) non-ME-associated piRNAs may have wider cellular roles. To explore this, predicted mRNA targets of differentially expressed piRNAs identified putative involvement in cancer-associated pathways. CONCLUSION: This study confirms previous miRNA observations, giving credence to our novel demonstration of global piRNA dysregulation in gonadal malignant GCTs, through both ME and non-ME-associated pathways, which likely contributes to GCT pathogenesis.

The in vivo RNA structurome of the malaria parasite Plasmodium falciparum, a protozoan with an A/U-rich transcriptome.

Plasmodium falciparum, a protozoan parasite and causative agent of human malaria, has one of the most A/T-biased genomes sequenced to date. This may give the genome and the transcriptome unusual structural features. Recent progress in sequencing techniques has made it possible to study the secondary structures of RNA molecules at the transcriptomic level. Thus, in this study we produced the in vivo RNA structurome of a protozoan parasite with a highly A/U-biased transcriptome. We showed that it is possible to probe the secondary structures of P. falciparum RNA molecules in vivo using two different chemical probes, and obtained structures for more than half of all transcripts in the transcriptome. These showed greater stability (lower free energy) than the same structures modelled in silico, and structural features appeared to influence translation efficiency and RNA decay. Finally, we compared the P. falciparum RNA structurome with the predicted RNA structurome of an A/U-balanced species, P. knowlesi, finding a bias towards lower overall transcript stability and more hairpins and multi-stem loops in P. falciparum. This unusual protozoan RNA structurome will provide a basis for similar studies in other protozoans and also in other unusual genomes.

Correction: MiR-277/4989 regulate transcriptional landscape during juvenile to adult transition in the parasitic helminth Schistosoma mansoni.

[This corrects the article DOI: 10.1371/journal.pntd.0005559.].

Smallpox vaccination induces a substantial increase in commensal skin bacteria that promote pathology and influence the host response.

Interactions between pathogens, host microbiota and the immune system influence many physiological and pathological processes. In the 20th century, widespread dermal vaccination with vaccinia virus (VACV) led to the eradication of smallpox but how VACV interacts with the microbiota and whether this influences the efficacy of vaccination are largely unknown. Here we report that intradermal vaccination with VACV induces a large increase in the number of commensal bacteria in infected tissue, which enhance recruitment of inflammatory cells, promote tissue damage and influence the host response. Treatment of vaccinated specific-pathogen-free (SPF) mice with antibiotic, or infection of genetically-matched germ-free (GF) animals caused smaller lesions without alteration in virus titre. Tissue damage correlated with enhanced neutrophil and T cell infiltration and levels of pro-inflammatory tissue cytokines and chemokines. One month after vaccination, GF and both groups of SPF mice had equal numbers of VACV-specific CD8+ T cells and were protected from disease induced by VACV challenge, despite lower levels of VACV-neutralising antibodies observed in GF animals. Thus, skin microbiota may provide an adjuvant-like stimulus during vaccination with VACV and influence the host response to vaccination.

RNA modifications detection by comparative Nanopore direct RNA sequencing.

RNA molecules undergo a vast array of chemical post-transcriptional modifications (PTMs) that can affect their structure and interaction properties. In recent years, a growing number of PTMs have been successfully mapped to the transcriptome using experimental approaches relying on high-throughput sequencing. Oxford Nanopore direct-RNA sequencing has been shown to be sensitive to RNA modifications. We developed and validated Nanocompore, a robust analytical framework that identifies modifications from these data. Our strategy compares an RNA sample of interest against a non-modified control sample, not requiring a training set and allowing the use of replicates. We show that Nanocompore can detect different RNA modifications with position accuracy in vitro, and we apply it to profile m6A in vivo in yeast and human RNAs, as well as in targeted non-coding RNAs. We confirm our results with orthogonal methods and provide novel insights on the co-occurrence of multiple modified residues on individual RNA molecules.

Short- and long-range cis interactions between integrated HPV genomes and cellular chromatin dysregulate host gene expression in early cervical carcinogenesis.

Development of cervical cancer is directly associated with integration of human papillomavirus (HPV) genomes into host chromosomes and subsequent modulation of HPV oncogene expression, which correlates with multi-layered epigenetic changes at the integrated HPV genomes. However, the process of integration itself and dysregulation of host gene expression at sites of integration in our model of HPV16 integrant clone natural selection has remained enigmatic. We now show, using a state-of-the-art 'HPV integrated site capture' (HISC) technique, that integration likely occurs through microhomology-mediated repair (MHMR) mechanisms via either a direct process, resulting in host sequence deletion (in our case, partially homozygously) or via a 'looping' mechanism by which flanking host regions become amplified. Furthermore, using our 'HPV16-specific Region Capture Hi-C' technique, we have determined that chromatin interactions between the integrated virus genome and host chromosomes, both at short- (<500 kbp) and long-range (>500 kbp), appear to drive local host gene dysregulation through the disruption of host:host interactions within (but not exceeding) host structures known as topologically associating domains (TADs). This mechanism of HPV-induced host gene expression modulation indicates that integration of virus genomes near to or within a 'cancer-causing gene' is not essential to influence their expression and that these modifications to genome interactions could have a major role in selection of HPV integrants at the early stage of cervical neoplastic progression.

Identification of Recurrent Mutations in the microRNA-Binding Sites of B-Cell Lymphoma-Associated Genes in Follicular Lymphoma.

Follicular lymphoma (FL) is a common indolent B-cell lymphoma that can transform into the more aggressive transformed FL (tFL). However, the molecular process driving this transformation is uncertain. In this work, we aimed to identify microRNA (miRNA)-binding sites recurrently mutated in follicular lymphoma patients, as well as in transformed FL patients. Using whole-genome sequencing data from FL tumors, we discovered 544 mutations located in bioinformatically predicted microRNA-binding sites. We then studied these specific regions using targeted sequencing in a cohort of 55 FL patients, found 16 recurrent mutations, and identified a further 69 variants. After filtering for QC, we identified 21 genes with mutated miRNA-binding sites that were also enriched for B-cell-associated genes by Gene Ontology. Over 40% of mutations identified in these genes were present exclusively in tFL patients. We validated the predicted miRNA-binding sites of five of the genes by luciferase assay and demonstrated that the identified mutations in BCL2 and EZH2 genes impaired the binding efficiency of miR-5008 and miR-144 and regulated the endogenous levels of messenger RNA (mRNA).

Locating patterns in Nanopore currents using time-warped signal representation of consensus nucleotides for demultiplexing and motif detection.

Nanopore-based approaches for the sequencing of DNA and RNA molecules are promising technologies with potential applications in clinical genomics. These approaches have generated large numbers of time series objects over the years, however, it remains a challenge to accurately decipher the underlying nucleotide sequence corresponding to a given signal. By using a combination of consensus signal averaging and stream monitoring of variable-length motifs, we outline an online pattern matching framework that can efficiently locate consensus sequences in real world Nanopore datasets. We demonstrate the applicability of our proposed framework across two use-cases: demultiplexing of DNA barcodes and multiple motif site identification in RNA transcripts.

3D model of harlequin ichthyosis reveals inflammatory therapeutic targets.

The biology of harlequin ichthyosis (HI), a devastating skin disorder caused by loss-of-function mutations in the gene ABCA12, is poorly understood, and to date, no satisfactory treatment has been developed. We sought to investigate pathomechanisms of HI that could lead to the identification of new treatments for improving patients' quality of life. In this study, RNA-Seq and functional assays were performed to define the effects of loss of ABCA12 using HI patient skin samples and an engineered CRISPR/Cas9 ABCA12 KO cell line. The HI living skin equivalent (3D model) recapitulated the HI skin phenotype. The cytokines IL-36α and IL-36γ were upregulated in HI skin, whereas the innate immune inhibitor IL-37 was strongly downregulated. We also identified STAT1 and its downstream target inducible nitric oxide synthase (NOS2) as being upregulated in the in vitro HI 3D model and HI patient skin samples. Inhibition of NOS2 using the inhibitor 1400W or the JAK inhibitor tofacitinib dramatically improved the in vitro HI phenotype by restoring the lipid barrier in the HI 3D model. Our study has identified dysregulated pathways in HI skin that are feasible therapeutic targets.

The midbody interactome reveals unexpected roles for PP1 phosphatases in cytokinesis.

The midbody is an organelle assembled at the intercellular bridge between the two daughter cells at the end of mitosis. It controls the final separation of the daughter cells and has been involved in cell fate, polarity, tissue organization, and cilium and lumen formation. Here, we report the characterization of the intricate midbody protein-protein interaction network (interactome), which identifies many previously unknown interactions and provides an extremely valuable resource for dissecting the multiple roles of the midbody. Initial analysis of this interactome revealed that PP1ß-MYPT1 phosphatase regulates microtubule dynamics in late cytokinesis and de-phosphorylates the kinesin component MKLP1/KIF23 of the centralspindlin complex. This de-phosphorylation antagonizes Aurora B kinase to modify the functions and interactions of centralspindlin in late cytokinesis. Our findings expand the repertoire of PP1 functions during mitosis and indicate that spatiotemporal changes in the distribution of kinases and counteracting phosphatases finely tune the activity of cytokinesis proteins.

Visualization and analysis of RNA-Seq assembly graphs.

RNA-Seq is a powerful transcriptome profiling technology enabling transcript discovery and quantification. Whilst most commonly used for gene-level quantification, the data can be used for the analysis of transcript isoforms. However, when the underlying transcript assemblies are complex, current visualization approaches can be limiting, with splicing events a challenge to interpret. Here, we report on the development of a graph-based visualization method as a complementary approach to understanding transcript diversity from short-read RNA-Seq data. Following the mapping of reads to a reference genome, a read-to-read comparison is performed on all reads mapping to a given gene, producing a weighted similarity matrix between reads. This is used to produce an RNA assembly graph, where nodes represent reads and edges similarity scores between them. The resulting graphs are visualized in 3D space to better appreciate their sometimes large and complex topology, with other information being overlaid on to nodes, e.g. transcript models. Here we demonstrate the utility of this approach, including the unusual structure of these graphs and how they can be used to identify issues in assembly, repetitive sequences within transcripts and splice variants. We believe this approach has the potential to significantly improve our understanding of transcript complexity.

Non-coding RNA Expression, Function, and Variation during Drosophila Embryogenesis.

Long non-coding RNAs (lncRNAs) can often function in the regulation of gene expression during development; however, their generality as essential regulators in developmental processes and organismal phenotypes remains unclear. Here, we performed a tailored investigation of lncRNA expression and function during Drosophila embryogenesis, interrogating multiple stages, tissue specificity, nuclear localization, and genetic backgrounds. Our results almost double the number of annotated lncRNAs expressed at these embryonic stages. lncRNA levels are generally positively correlated with those of their neighboring genes, with little evidence of transcriptional interference. Using fluorescent in situ hybridization, we report the spatiotemporal expression of 15 new lncRNAs, revealing very dynamic tissue-specific patterns. Despite this, deletion of selected lncRNA genes had no obvious developmental defects or effects on viability under standard and stressed conditions. However, two lncRNA deletions resulted in modest expression changes of a small number of genes, suggesting that they fine-tune expression of non-essential genes. Several lncRNAs have strain-specific expression, indicating that they are not fixed within the population. This intra-species variation across genetic backgrounds may thereby be a useful tool to distinguish rapidly evolving lncRNAs with as yet non-essential roles.

Terminal uridylyltransferases target RNA viruses as part of the innate immune system.

RNA viruses are a major threat to animals and plants. RNA interference (RNAi) and the interferon response provide innate antiviral defense against RNA viruses. Here, we performed a large-scale screen using Caenorhabditis elegans and its natural pathogen the Orsay virus (OrV), and we identified cde-1 as important for antiviral defense. CDE-1 is a homolog of the mammalian TUT4 and TUT7 terminal uridylyltransferases (collectively called TUT4(7)); its catalytic activity is required for its antiviral function. CDE-1 uridylates the 3' end of the OrV RNA genome and promotes its degradation in a manner independent of the RNAi pathway. Likewise, TUT4(7) enzymes uridylate influenza A virus (IAV) mRNAs in mammalian cells. Deletion of TUT4(7) leads to increased IAV mRNA and protein levels. Collectively, these data implicate 3'-terminal uridylation of viral RNAs as a conserved antiviral defense mechanism.

Defective germline reprogramming rewires the spermatogonial transcriptome.

Defective germline reprogramming in Piwil4 (Miwi2)- and Dnmt3l-deficient mice results in the failure to reestablish transposon silencing, meiotic arrest and progressive loss of spermatogonia. Here we sought to understand the molecular basis for this spermatogonial dysfunction. Through a combination of imaging, conditional genetics and transcriptome analysis, we demonstrate that germ cell elimination in the respective mutants arises as a result of defective de novo genome methylation during reprogramming rather than because of a function for the respective factors within spermatogonia. In both Miwi2-/- and Dnmt3l-/- spermatogonia, the intracisternal-A particle (IAP) family of endogenous retroviruses is derepressed, but, in contrast to meiotic cells, DNA damage is not observed. Instead, we find that unmethylated IAP promoters rewire the spermatogonial transcriptome by driving expression of neighboring genes. Finally, spermatogonial numbers, proliferation and differentiation are altered in Miwi2-/- and Dnmt3l-/- mice. In summary, defective reprogramming deregulates the spermatogonial transcriptome and may underlie spermatogonial dysfunction.

Genomic positional conservation identifies topological anchor point RNAs linked to developmental loci.

BACKGROUND: The mammalian genome is transcribed into large numbers of long noncoding RNAs (lncRNAs), but the definition of functional lncRNA groups has proven difficult, partly due to their low sequence conservation and lack of identified shared properties. Here we consider promoter conservation and positional conservation as indicators of functional commonality. RESULTS: We identify 665 conserved lncRNA promoters in mouse and human that are preserved in genomic position relative to orthologous coding genes. These positionally conserved lncRNA genes are primarily associated with developmental transcription factor loci with which they are coexpressed in a tissue-specific manner. Over half of positionally conserved RNAs in this set are linked to chromatin organization structures, overlapping binding sites for the CTCF chromatin organiser and located at chromatin loop anchor points and borders of topologically associating domains (TADs). We define these RNAs as topological anchor point RNAs (tapRNAs). Characterization of these noncoding RNAs and their associated coding genes shows that they are functionally connected: they regulate each other's expression and influence the metastatic phenotype of cancer cells in vitro in a similar fashion. Furthermore, we find that tapRNAs contain conserved sequence domains that are enriched in motifs for zinc finger domain-containing RNA-binding proteins and transcription factors, whose binding sites are found mutated in cancers. CONCLUSIONS: This work leverages positional conservation to identify lncRNAs with potential importance in genome organization, development and disease. The evidence that many developmental transcription factors are physically and functionally connected to lncRNAs represents an exciting stepping-stone to further our understanding of genome regulation.

MicroRNA degradation by a conserved target RNA regulates animal behavior.

microRNAs (miRNAs) repress target transcripts through partial complementarity. By contrast, highly complementary miRNA-binding sites within viral and artificially engineered transcripts induce miRNA degradation in vitro and in cell lines. Here, we show that a genome-encoded transcript harboring a near-perfect and deeply conserved miRNA-binding site for miR-29 controls zebrafish and mouse behavior. This transcript originated in basal vertebrates as a long noncoding RNA (lncRNA) and evolved to the protein-coding gene NREP in mammals, where the miR-29-binding site is located within the 3' UTR. We show that the near-perfect miRNA site selectively triggers miR-29b destabilization through 3' trimming and restricts its spatial expression in the cerebellum. Genetic disruption of the miR-29 site within mouse Nrep results in ectopic expression of cerebellar miR-29b and impaired coordination and motor learning. Thus, we demonstrate an endogenous target-RNA-directed miRNA degradation event and its requirement for animal behavior.